Spectral Properties of the Two-Dimensional Laplacian with a Finite Number of Point 

Interactions 



O 

o 
o 



T. Shigehara^, H. Mizoguchi-^, T. Mishima^ and Taksu Cheon^ 
^Department of Information and Computer Sciences, Saitama University 
Shimo-Okubo, Urawa, Saitama 338 Japan 
^Laboratory of Physics, Kochi University of Technology 
Tosa Yamada, Kochi 782 Japan 

We discuss spectral properties of the Laplacian with multiple (A'^) point interactions in two- 
dimensional bounded regions. A mathematically sound formulation for the problem is given within 
the framework of the self-adjoint extension of a symmetric (Hermitian) operator in functional anal- 
ysis. The eigenvalues of this system are obtained as the poles of a transition matrix which has size 
■ A*'. Closely examining a generic behavior of the eigenvalues of the transition matrix as a function of 

, the energy, we deduce the general condition under which point interactions have a substantial effect 

' on statistical properties of the spectrum. 

^ ! 
> ■ 

O ■ I. INTRODUCTION 

O: 

I I The Laplacian is one of fundamental operators which describe real systems and indeed appears in various fields 
such as classical mechanics, electromagnetic dynamics, fluid dynamics, quantum mechanics and so on. The eigenvalue 
y—i ' problem of the Laplacian, in particular, has a direct relation to the microscopic world where quantum mechanics 
^ : governs the dynamics. The eigenvalues of the Laplacian on a bounded region, uji, i = 1, • • •, exactly correspond to 
1^ ' the allowable energies of a free particle moving in a microscopic billiard, which are discrete in general, and each 
eigenf unction, il>i{x), describes the wave function of the corresponding eigenstate, the square of which, |'0i(a?)P, 
expresses the probability such that the particle with energy uji exists at a position x in the billiard. 

The quantum billiard problem is a natural idealization of the particle motion in microscopic bounded regions. 
[ — I The one-electron problem in "quantum dots" is a possible setting which may be used as a single-electron memory, a 
0^ ' promising computational device in the future. It is now possible to actually construct such structures with extremely 
, pure semiconductors thanks to the rapid progress in the mesoscopic technology. Real systems are, however, not free 
^ ^' from impurities which affect the particle motion inside. In the presence of a small amount of contamination, even a 
^ I ^ ' single-particle problem becomes unmanageable. The modeling of the impurities with point interactions is expected 
to make the problem easy to handle without changing essential dynamics. 

Apart from the applicability to physical systems with a microscopic or mesoscopic scale, the Laplacian with point 
^ , interactions is an attractive object from a viewpoint of dynamical systems. The two-dimensional quantum billiard is an 
appropriate tool for examining generic features of dynamical systems because of the wide rage of dynamical behaviors, 
going from the most regular (integrable) to the most irregular (chaotic) depending on the geometry of its boundary. 
Although no mathematical proof exists, it is widely believed that fingerprints of the regular or irregular nature of 
^ . , the classical motion can be found in statistical properties of the spectrum (eigenvalues and eigenfunctions) in the 
' corresponding quantum system |^ . One of such statistical measures is the nearest-neighbor level spacing distribution 
P{S), which is defined such that P{S)dS is the probability to find the spacing between any two neighboring energy 
levels in the interval (5, S + dS). Integrable systems such as circular, elliptic and rectangular billiards obey Poisson 
statistics; 

P{S)^cxpi-S). (1) 

(It is conventional to take a unit of the energy such that the average level spacing is one. Namely, P{S) satisfies 
SP{S)dS = 1 as well as P{S)dS = 1.) On the other hand, chaotic systems such as Sinai's billiard and 
Bunimovich's stadium are described by the prediction of the Gaussian orthogonal ensembles; 

ttS ( 7r5^\ 

P(5) = — exp -— . (2) 
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There are other statistics suitable for a measure of the degree of regularity (or irregularity), such as As-statistic 
introduced by Dyson and Mehta [^,^. From a viewpoint of dynamical systems, the integrable billiards with point 
interactions inside are classified into a category called "pseudointegrable" . The nature of classical motion in the 
pseudointegrable systems is integrable in the sense that a set of unstable trajectories are of measure zero in the phase 
space. However, several numerical experiments show that under a certain condition, quantization induces the chaotic 
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energy spectra, which might be regarded as a counterexample for the correspondence between the level statistics and 
the underlying classical motion. 

One of the main purposes of this paper is to clarify the condition under which the statistical properties of quan- 
tum spectrum are substantially affected by point interactions. We spend a large portion of the paper in giving a 
mathematical background to answer this question. In spite of its seeming simplicity, a careful treatment is demanded 
for point interactions in quantum mechanics in case of the spatial dimension d > 2. Relying on the self-adjoint 
extension theory in functional analysis, we give a general formula suitable for the resolvent (Green's function) which 
describes a particle propagation in quantum billiard with point interactions. It is hard to proceed further without 
losing any generality; Numerics is, in general, demanded to calculate each eigenvalue (pole of the Green's function) 
and the corresponding eigenfunction. By introducing suitable approximations for examining statistical properties on 
spectrum, however, we deduce the condition for the strong coupling where the effect of point interactions on the 
quantum spectrum is maximal. 

The paper is organized as follows. In Sec. II, the Green's function for the Laplacian with a finite number of point 
interactions is deduced in a mathematically sound manner. We discuss the condition for the strong coupling for the 
case of a single interaction in Sec. Ill, which is extended to the case of multiple interactions in Sec. IV. The present 
work is summarized in Sec.V. 



II. FORMULATION OF QUANTUM BILLIARDS WITH POINT INTERACTIONS 

We start from an empty billiard. Let us consider a quantum particle of mass M moving freely in a bounded region 
S in two spatial dimension. The wave functions are assumed to vanish on its boundary dS. We denote the eigenvalues 
and corresponding eigenfunctions of the system as and ipn, namely 

i^O = Hoifnix) = Snfnix), n = 1, • • • (3) 

with 

{x) = where x G dS. (4) 

The Hamiltonian Hq is the kinetic operator in LP'{S) with domain D{Ho) = H^{S) fl Hq{S) in terms of the Sobolev 
spaces. According to the Weyl formula, the average level (eigenvalue) density is given by 

Pa. = (5) 

where we denote the area of S by the same symbol. Note that pav is independent of the energy w in two dimension. 
Assuming (pn {x) to be normalized to unity, the Green's function of Hq is given by 

G(0)(f,f ;a;) = (a; - Ho)-\S,^) = V '^"(^)'^"(^') . (6) 

„=i ~ ^" 

We now place N point interactions at xi,...,xn- Naively, one defines the interactions in terms of the Dirac's delta 
function in two dimension; 

JV 

i=l 

However, the Hamiltonian H is not mathematically sound for spatial dimension d>2. This can be seen even in case 
of a single interaction (iV = 1). In this case, the eigenvalue equation of H is reduced to 

V ^En^ = v-\ (8) 

However, since the average level density is constant with respect to the energy, the infinite series on the left hand 
side does not converge in two dimension. For higher dimensions, the average level density has a positive power 
dependence of the energy, which makes the series divergent. To handle the divergence, a scheme for regularization 
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and renormalization is called for. One of the most mathematically satisfying schemes is given by the self-adjoint 
extension theory of functional analysis [Q . 

We first consider in L^(S') the nonnegative operator 



(9) 

C-(5-X) 



with its closure Hx in L^(S'), where we set X — {xi, ...,xn} ■ Namely, we restrict D{Hq) to the functions which 
vanish at the position of the point interactions. By using integration by parts, it is easy to prove that the operator 
Hx is symmetric (Hermitian). But it is not self-adjoint. Indeed, the equation 

H^^P{lu,x)^uj'iP{uj,x), %l^eD{H*x), Imu^O, (10) 

has the N independent solutions 

ip,{Lu,x) ^ G'-°\x,x,;uj), xeS-X, i = l,---,iV, (11) 

indicating 



D{H*x) = D{Hx) ® Ran{Hx - uj) ® Ran{Hx - u)) ^ D{Hx), Imuj^Q. (12) 

Since Hx has the deficiency indices {N,N), Hx has, in general, A^^-parameter family of self-adjoint extensions. All 
self-adjoint extensions Hu^x of Hx are given by 



N N 

D{Hu,x) = {f + Y. ^»(^^+ + E U,ji;3-)\f ^Hx,c,e C}, 
i=i j=i 

N N _ N N 

Hu,x{f + E + E U^3^J-)} = Hxf + zA ^ C.(V.+ " E t^yV'j-), (13) 

1=1 3 = 1 i=l j=l 

where Uij denotes a A^-dimensional unitary matrix and 

^j±{x)^^,{±tA,x), A>0, j = l,---,7V, (14) 



provide a basis for Ker{H^ =F iK) = Ran{Hx ± *A) , respectively. Positive A is regarded as a mass scale which can 
be arbitrarily fixed. The operator Hu^x corresponds to the Hamiltonian for the system with N point interactions. 
The special case U — —1 leads to the kinetic operator in L'^{S), H^i^x — Hq, since — ipj^ E D{Ho). 

By using Krein's formula, we can obtain the relation between two resolvents {uj — Hq)~^ and {to — Hjjx)^^ for the 
unperturbed and perturbed systems respectively. Assume that v4 is a densely defined, closed symmetric operator in 
some Hilbert space with deficiency indices {N, N). Let B and C be two self-adjoint extensions of A and denote by A 
the maximal common part of B and C. Let M, < M < N, be the deficiency indices of yl and let {(pi{Lu), • • • , (/)7\/(a;)}, 
which are linearly independent, span the corresponding deficiency subspace of A; 

A*(t>,{uj) ^ Lj(t).,{Lj), (t).,eD{A*), i = l,---,M, Imuj^O. (15) 

Then Krein's Formula reads 

Theorem Let B, C, A and A be as above. Then 

M 

(^-B)-i-(a>-C)-^ =- ^ A,,(c^)(</),(cD),.)</),;(c^), u;ep{B)np{C), (16) 

where the matrix X{lu) is nonsingular for u € p{B) n p{C) and Ay (w) and 0j(aj), i, j — 1, • ■ • , M, may be chosen to 
be analytic in a; G p{B) n p{C). Indeed, 4>i{Lo) can be defined as 

(P,{uJ)=(j)^{uJo)-{uJ-UJo){uJ-C)-^(p^{uJo), i^l,---,M, LJ E p{C) , (17) 



where (/)i{ujo), i — 1, • • • ,M, Im cjq 7^ 0, are linearly independent solutions of Eq.(|l^) for uj — uo^) and the inverse 
matrix of A(ijj) satisfies 
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X-^{cj) = X-^{u;')-{u;-u;'){cl>,{Q),M^')}, i,j = l,---,M, cu , u;' € p{B) n p{C) , (18) 
if the (j)i{uj), i = 1, ■ ■ ■ , M , are defined according to Eq.(|T] 



Applying Krein's formula {A — Hx, B — Hjj x, C — Hq and hence A — Hx, M — N), we obtain the relation for 

N 

(w - Hu,x)-^ = (co - Hor^ - ^^3i^){i'ji^)^ ■)M^), (19) 

where A~^(a;) satisfies 

A-.i(c.) - A-.i(^') = -{lu - c.')(^,(t^),^.(^')> = Gg^(^) - Gfiu;'), z,j = 1, • • • ,iV. (20) 

Here we set 

Gf(u.) ^ G(°)(x.,.-;c.) = f (21) 

The second equality in Eq.(po|) follows from the resolvent formula 

{lu ~ Lu'){iu - Ho)-\u;' - Ho)-' = {uj' - Ho)-' - (lu - Ho)-\ (22) 

Each term on the right hand side in Eq.(po|) is a divergent series for i = j. However, the divergence cancels each 
other. The most important indication of Krein's formula is that it is sufficient to define A(a-') for some fixed lu since 
then A(a;) for any uj follows from Eq.(pO|). 

The Green's function of Hu^x is given by the coordinate representation of the resolvent Eq.(p^); 

N 

Gu^x{x,x';u;) = {lo - Hu,x)'HS, x') ^ G^°\x, x' -,0;) - ^ G^°\x,xf,u;)K,{u;)G^°\x,,x';u;). (23) 

From Eq.(|2^), we recognize that — A(w) corresponds to the transition matrix (T-matrix) in the presence of N point 
interactions. It follows from Eq.(n3[) that 



AT N 

{oj - Hu^x){^i+ + Uij^pj^) = {uj- iA)i^,+ + {uj + iA) J2 Uijtpj-. (24) 
Setting u) — — iA in Eq.(p^), we get 

1 ^ 

(-zA - Hu,x)-'i^^+ = -^(^«+ + E U^i^J-)- (25) 
On the other hand. Eq.([l9|) shows that the left hand side in Eq.([25|) is written as 

N 

{~iA- Hu,x)-'A+ - HA--Ho)-'A+ ~ E XJk{-^A){^k+,i^^+H3-. (26) 

3:k=l 

The resolvent formula (E2h indicates 



(-zA - Ho)-'A+ = -^(V'.+ - ^^-), (27) 

= J^(Gi';)(-zA) - G^m) = ±.iX-li-^A) - A^/(zA)), (28) 

where the second equality in Eq.(p^) follows from Eq.([20|). Substituting Eqs.(p7|) and ( p8| ) into Eq.(|2^) and compared 
with Eq.(|2q), we obtain the relation between the unitary matrix and the T-matrix, 
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N N 

U.J = -6,j + J2 A,fe(-*A)(A,-/(-*A) - A,-/(*A)) = - ^ A,,(-»A)A^/(zA), (29) 

k=l k=l 

or equivalently 

[/ = -*[A(-iA)A-i(iA)] (30) 
in a simple matrix form. Noticing that Eq.([ig|) implies 

A(a;)t = X{oj), (31) 
we realize that the unitarity of U is equivalent to the fact that A(?A) (resp. A(— iA)) is a normal matrix. 

III. THE CASE OF SINGLE INTERACTION 

In case of a single point interaction, a general form of X{uj) G C which satisfies Eqs.(|20|) and (|3^) is given by 

X-\uj) ^ G{uj) - v-\ (32) 

where 

G(c.) = ^«(-^i)'(^73- + 72Ta^)' (33) 
.-=Acot^5:gil^, 0<^<2.. (34) 



(Actually it would have been sufficient to define A(iA) since Xjuj) is calculated from Eq.(gO|) with = 1, as mentioned 
in the previous section.) At first sight, the second term in Eq.(|33[) seems to be redundant since it is energy-independent. 
However, it plays an essential role in making the problem well-defined; The infinite series in Eq.(|33|) does not converge 
without the second term. For spatial dimension d > 4, in general, the summation in Eq.(|3^) diverges. This reflects 
the fact that the billiard problem with point interactions is not well-defined for d > 4 in quantum mechanics. The 



energy independence of the counterterm also ensures that A(w) in Eq.(|32|) satisfies Eq.(20). (Note that, apart from the 

counterterm, G{uj) has the same form as G^ii{uj).) As one varies < 9 < 2it for some fixed A, the value of v, which 
is also energy-independent, ranges over the whole real number. We can formally consider z; as a coupling strength of 
the point interaction. Its relation to physical strength of the interaction has been discussed elsewhere m. Here we 
just mention that the system approaches the empty billiard in the limit oi v ^ 0. Inserting A(ix') in Eq.(p3) together 
with Eqs.(|3^) and ( ^4[ ) into Eq.(^), we get a simple expression of U which manifests the unitarity; 

t/ = -e*^ (35) 

irrespective to the value of A. 

The eigenvalues for the system perturbed by a single point interaction with coupling strength v located at xi are 
obtained as the solutions of the equation 

X-\uj) = 0, (36) 

namely 

G{uj)=v-\ (37) 

Eq.(p3[) shows that within any interval between two neighboring unperturbed eigenvalues, G{uj) is a monotonically 
decreasing function that ranges over the whole real number. This means that Eq.(^^ has a single solution on each 
interval; 

El < uji < e2 < UJ2 < £3 < ^3 < , for V > 0, , . 

uii < El < UJ2 < £2 < ^^3 < £3 < , for w < 0. ^ ' 



5 



The eigenfunction corresponding to an eigenvalue aj„ is written in terms of Green's function of the empty biUiard as 



k=l 



— £k 



(39) 



with a normahzation constant 




(40) 



In order to obtain each solution of Eq.(|3^), a numerical task is needed in general. Our main purpose is, however, 
to examine under what conditions a point interaction brings about a significant effect in the statistical properties of 
spectrum. To achieve our aim without losing the essence, while still keeping loss of generality minimal, we introduce 
some approximations as follows. The first (and main) one is that the value of (the square of) each wave function at 
xi is replaced by its average among many; 



(p„(fi)^ ~ {(pnixif) = 1/S. 



(41) 



Since it is often the case that the statistics are taken within a large number of, sometimes thousands of eigenstates 
(with xi and v fixed), Eq.(^l|) is quite satisfactory. Note that the second equality in Eq.(^) is valid irrespective to 
the exact position of the interaction xi as well as the energy region where the average is taken. Keeping Eq.(|4l|) 
in mind, we recognize from Eq.(^9|) that if a perturbed eigenvalue a;„ is close to an unperturbed one e„ (or £„+i)j 
then the corresponding eigenfunction ipn is not substantially different from cp„ (or Lpn+i)- Thus, the disturbance 
by a point interaction is mainly restricted to eigenstates with an eigenvalue around which G(ct;) has an inflection 
point. This is because each inflection point of G{lu) is expected to appear, on average, around the midpoint on the 
interval between two neighboring unperturbed eigenvalues. Let (a),G(ii')) be one of such inflection points of G{lu); 
uj ~ (e,„ + £m+i)/2 for some m. In this case, the contributions on G{Cj) from the terms with n m cancel each 
other, and we can approximate the summation in Eq.(p3|) by a principal integral with a high degree of accuracy; 



Giiv)^9iiv), 



1 



E 



-E E'^ + K^ 



dE, 



(42) 
(43) 



where we have defined a continuous function g(uj) which behaves like an interpolation of the inflection points of G{lu). 
Clearly, G{uj) ~ ^(w) is valid only aromid the inflection points of G{uj). Using an elementary indefinite integral 



1 



E 



E E^ +A'^ 



dE = - In ■ 



\^-E\ 



we obtain 



G{UJ) ~ {LPn{Xi) )pai.ln^ ~ — In - , 



(44) 



(45) 



where the second equality follows from Eqs.(H) and (^l|). Eq.(^) indicates that the maximal coupling of a point 
interaction is attained with the coupling strength v which satisfies 



M , oj 

t:— In -r- 

2tt a 



(46) 



The "width" of the strong coupling region (allowable error of w ^ in Eg . mm ) can be estimated by considering a 
linearized eigenvalue equation. Expanding G(w) at w = a), we can rewrite EqT(p7|) as 



G{Cj) + G' {Cj){uj - Cj) ~ V- 



(47) 



G'(w)(a 



M , w 

t:— In -r- 

27r A 



(48) 
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In order to ensure that the perturbed eigenvalue is close to w, the range of the right hand side in Eq. 
restricted to 



M , bj 

7;— In -r 
27r A 



< 



has to be 



(49) 



where the width 6^-1 is defined by 



S^-i = |G'(J))| pj. 



(50) 



This is nothing but the average variance of the linearized G on the interval between the two unperturbed eigenvalues 
just below and above Cj. Using the approximation in Eq.(B^), the value of I G" ((D) I can be estimated as follows; 



= 8((p„(xi)2)p2^£ 



1 



1 



2 

av ' 



(51) 



The second equality follows from the approximation that the unperturbed eigenvalues are distributed with a mean 
interval p~y in the whole energy region. This assumption is quite satisfactory, since the denominator of G'{llj) is of 



the order of (w — e„ 



indicating that the summation in Eq.(pT|) converges rapidly. From Eq.(|T|), we obtain 

^c-i ~ 7r2((^„(fi)2)p^„ ~ (52) 



From Eqs.(|49|) and (p^), we recognize that the effect of a point interaction of coupling strength v on statistical 
properties of spectrum is substantial only in the eigenstates with eigenvalue uj such that 



M , UJ 
— In — 
27r A 



< 



(5,1-1 



M. 



(53) 



We do not go into details of numerical experiments in this paper, but just mention that the validity of the conjecture 
( p3| ) has been confirmed by examining the spectrum of rectangular billiards with a single point interaction inside 
Concerning the nearest-neighbor level spacing distribution P{S), a level repulsion seen in Eq.(|^) is observed under 
the condition (p3|), while P{S) is not substantially different from the Poisson prediction, Eq.(|l|), otherwise. 

Before closing this section, we give a short comment on the shape of the billiard. Our implicit assumption for the 
shape is that the average level density of the empty billiard is dominated by the area term, which does not depend 
on energy. The assumption is justified for a generic two-dimensional billiard which has the same order of length scale 
in each direction, irrespective to a full detail of the shape of the billiard. 



IV. THE CASE OF MULTIPLE INTERACTIONS 

We intend to generalize Eq.(|3^) to finitely many point interactions. This is achieved by introducing the iV-parameter 
family of self-adjoint extensions of Hx defined by the unitary matrix U, Eq. (|30|) , together with 

,-1, X j Gi{uj)~v-\ for i=j, VieR, . 

^-(")=\G(f(c.), forz^,, ^,J = lr--,iV, (54) 

where G^^\uj) is defined in Eq.(^l[) and 

00 ^ 

G.H = E^"(^')'(— — + ;:jTa^)' ^ = h---,N. (55) 

Eq.(p4[) with opposite sign describes the T-matrix of the system with N point interactions which satisfy separated 
boundary conditions at each point Xi . The coupling strength of the i-th interaction is assigned by a real number Vi , 
I = 1, • • • , A''. The eigenvalues of this system are the poles of the T-matrix, determined by 

det X'\lo) = 0. (56) 
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Qualitative behavior of the eigenvahics of A ^{u;) as a function of w can be examined as follows. Let us cut the infinite 
series of each matrix element of A^^(w) at Umax which satisfies uj <^ ^umax- then have 



X-\uj) = lim 

rimax^OO 

where is a sum of an "unperturbed part" T^^^ and an "interaction part" V^^\io) 

E(0)(a;)=r(0) + y(°)(w) 



(57) 



(58) 



with 



J^(O) ^ 



Here we set 







AO) 







n=l 



(^f(fi) ••• ^^(fiv))- (59) 



"'max 



n=l 



A2 



(60) 



The iV-dimensional matrix Y^-^^iuS) formally has a form similar to the Hamiltonian with Umax rank-one (separable) 
interactions. We realize that takes a role of the "unperturbed energy" , while the inverse of — corresponds 
to the "strength" of the n-th interaction. In the following, we assume that < d'^''^ < • • • < rf^' without a loss of 
generality. 

The eigenvalues of S(°)(a;) can be obtained in a perturbative manner as follows. We start by setting k = 1 and 
define 



UJ-Sk 



(61) 



Namely, only the first term in the summation of V^'^ '^\'^) is taken into account in the matrix T'-'^\lij). Since T'^^^uj) 

has a single rank-one interaction, we obtain the eigenvalues of T^'^\ijj:) as the solutions d^^\ i = 1,---,N, of the 
equation 



N (fe-1) 



E 



(62) 



It is worthy to note that 

dfi'' < < d^-'^ for u; < e^, rff < rff ^ < ^^^'^ for cj > Sk. (63) 
We define an orthogonal transformation Q^''^ (w) such that the matrix T^*^) (w) is diagonal in the new basis. 



/d['^ 



T('=)(a;) =*n('=)(a;)T('=)(a;)0('=)(a;) = 











(64) 



Clearly, the i-th column of ^^''^uj) is given by the normalized cigcnfunction corresponding to the z-th eigenvalue d^^^\ 
the exact form of which is unnecessary in the present argument. With the new basis, the new Hamiltonian matrix 
has a form 
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where 



(65) 



n=i;+l 



UJ - Er 



(66) 



together with 

(^i')(xi) ... ^l'=)(fA.)) = (^i'^-'^(xi) ... ^ri)(x^))f7W(c.), 



7 ""max- 



(67) 



The new interaction matrix y (oj) keeps the separabihty of the interaction, while the number of the terms decreases 
by one, compared to V''^~^\uj). The procedure mentioned above can be repeated for k = 2, - ■ ■ ,nmax, successively. 
After the repetition, we obtain i = 1, ■ ■ ■ , N , which are exactly the eigenvalues of E'^°)(w). 

The above algorithm serves to understand the qualitative behavior of the eigenvalues of A~^(a;). The eigenvalue 
equation (p3) indicates that if Sk is far from the energy oj under consideration, the "strength" of the rank-one 
interaction7(a; — £k)~^, is very weak and it has little effect except that each eigenvalue shifts by 



(fc-l) / ^ \2 

/fc) _ ^(fc-1) _^ fk y^i) _ ^(fc-i) 

to — Sk 



(68) 



This means that fl^'''>{u;) substantially differs from the unit matrix only in case of Sk — Thus, we can approximate 
l](0)(w) by 



E(")(w)~r(w) + ^(a;) 



(69) 



with 







Here we set 



/d. \ 

r i m a. a: 



( ipn{xi) ■ ■ ■ (Pn{xN) ) 



(70) 



•A2 



E 

n=l 
En 



^n{XiY 



uj - e„ 



z = l,...,7V. 



(71) 



Note that di shows the tracks of the discard terms in ^(a;). In the hmit of 

di 



oo, we can estimate 



— In V- 

27r A 



(72) 



as before. Since the shift of eigenvalues caused by rank-one interactions is limited [see Eq.(|6^)], we expect from 
Eq.([70|) together with Eq.(^) that the condition ( p3| ) for a single point interaction can be generalized to the case of 
multiple interactions. Namely, the i-th interaction with strength u,; affects the spectral properties in the energy region 
which satisfies 



V. CONCLUSION 



We have discussed the spectral properties ol two-dimensional Laplacian with a finite number of point interactions. 
In spite of the apparent simphcity, careful treatments are required for point interactions. Based on the self-adjoint 
extension theory in functional analysis, we have deduced the Green's function appropriate for the operator. The 
problem is based on obvious physical motivations. The operator exactly corresponds to the quantum-mechanical 
Hamiltonian for a point particle moving in a two-dimensional bounded region (billiard) with point impurities inside. 
Based on a general argument, we have clarified the condition under which the statistical properties of spectrum are 
substantially influenced by the point interactions. The findings are summarized as follows; 

(1) For a two-dimensional billiard, the effect of a point interaction with coupling strength v on statistical properties 
of spectrum is maximal under the condition 

where M is the mass of a particle moving in the billiard and A is an arbitrary mass scale. This indicates that the 
maximal coupling region shifts with a logarithmic dependence of the energy oj in two dimension. 

(2) The width (J^-i (or an allowable error in v'"^ to look for the effect) is estimated as 



(75) 



which is energy- independent. If the value oi v differs from the right hand side in Eq.(^4|) to the extent of Eq.([75|) 
at energy uj, the effect of the point interaction on the statistics tends to disappear. 

(3) (1) and (2) are generalized to the case of multiple point interactions under separated boundary conditions. If we 
collect the interactions with the same order of magnitude of the coupling strength as a single group, we expect that 
the interactions belonging to one of such groups disturb the particle motion in a "coherent" manner in the energy 
region determined by Eq. ( fz^ , while their influence hardly appears in different energy regions where the difference 
between both sides in Eq.(|7^) is larger than Eq. ([zsj) . 
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